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!>- ! Abstract 

We present a search for the decay B + — > (£ = r, /i, or e) in (458.9±5.1) x 10 6 T(4S) decays recorded 
with the BXBAR detector at the SLAC PEP-II B -Factory. A sample of events with one reconstructed exclusive 
O ■ semi-leptonic B decay (B~ — > D°£~z/X) is selected, and in the recoil a search for B + — > £ + vi signal is 
performed. The r is identified in the following channels: r + — > e + v e V T , r + — > v^p T , r + — > n + V T , 
and r + — > ir + ir T' T . The analysis strategy and the statistical procedure is set up for branching fraction 
extraction or upper limit determination. We determine from the dataset a preliminary measurement of 
B(B + -> t + Vt) = (1-8 ± 0.8 ± 0.1) x 10~ 4 , which excludes zero at 2Ao, and f B = 230 ± 57MeV. 
Combination with the hadronically tagged measurement yields B(B + — > t + u t ) = (1.8 ± 0.6) x 10~ 4 . 
We also set preliminary limits on the branching fractions at B{B + — > e + v e ) < 7.7 x 10 -6 (90% C.L.), 
B(B + -> n+Vn) < 11 x 10~ 6 (90% C.L.), and B(B + -> r+^ r ) < 3.2 x 10" 4 (90% C.L.). 
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1 INTRODUCTION 



In the Standard Model (SM), the purely leptonic decay B + — ► £ + U£ 8 proceeds via quark annihilation into 
a W + boson (Fig. 1). 

The branching fraction is given by: 

2 

WllK fe | 2 , (i) 

where we have set Ti = c = 1, Gp is the Fermi constant, V u b is a quark mixing matrix element [1, 2], 
fs is the B + meson decay constant, which describes the overlap of the quark wave-functions inside the 
meson, t b + is the B + lifetime, and m B and m T are the B + meson and r masses. This expression is entirely 
analogous to that for pion decay. Physics beyond the SM, such as two-Higgs doublet models, could enhance 
or suppress the B(B + — ► t + ve) through the introduction of a charged Higgs boson [3]. 

Current theoretical values for f B (obtained from lattice QCD calculations) [4] have large uncertainties, 
and purely leptonic decays of the B + meson may be the only clean experimental method of measuring f B 
precisely. Given measurements of |V^{,| from semi-leptonic B — > ulv decays, f B could be extracted from 
the measurement of the B + — ► t + v t branching fraction. In addition, by combining the branching fraction 
measurement with results from B mixing, the ratio |V^&|/| V t d\ can be extracted from B{B + — > t + v t ) / Am, 
where Am is the mass difference between the heavy and light neutral B meson states. 

b vi 




u £ + 

Figure 1: The purely leptonic B decay B + — ► t + v t proceeding via quark annihilation into a W + boson. 

The SM estimates of these branching fraction are B(B + -> t + v t ) = (1.2 ± 0.4) x 10' 4 , B{B + -> 
= (5.6 ± 1.7) x 10^ 7 , and B(B + e+i/ e ) = (1.3 ± 0.4) x 10 -11 . We use \V ub \ = (4.43 ± 
0.54) x 10 -3 [5] and a theoretical calculation of / B = 189 ± 27 MeV [4] in Eq. 1. The differences in the 
branching fractions are due to helicity supression, which is expressed in Equation 1 via the different masses 
of the leptons. 

8 Charge-conjugate modes are implied throughout this paper. The signal B will be denoted as a B + decay while the semi- 
leptonic B will be denoted asaB". 



G\m B m\ 
8^ 



ml 



8 



2 THE BaBAR DETECTOR AND DATASET 



The data used in this analysis were collected with the BaBAR detector at the PEP-II storage ring. The sample 
corresponds to an integrated luminosity of 417.6 fb _1 with center-of-mass energy equal to the T(4S) rest 
mass (on-resonance) and 42.2 fb -1 taken 40MeV below BB threshold (off -resonance). The on-resonance 
sample consists of about (458.9 ± 5.1) x 10 6 T(4S) decays (BB pairs). The collider is operated with 
asymmetric beam energies, producing a boost of f3~j « 0.56 of the T(4S) along the collision axis. 

The BaBAR detector is optimized for asymmetric energy collisions at a center-of-mass (CM) energy 
corresponding to the T(4S) resonance. The detector is described in detail in Ref. [6]. The components 
used in this analysis are the tracking system composed of a five-layer silicon vertex detector and a 40-layer 
drift chamber, the Cherenkov detector for charged ir-K discrimination, the Csl calorimeter for photon and 
electron identification, and the flux return located outside of the 1.5T solenoidal coil and instrumented with 
resistive plate chambers for muon and neutral hadron identification. For the most recent 121 fb _1 of data, a 
portion of the muon system has been upgraded to limited streamer tubes [7, 8]. 

A GEANT4-based [9] Monte Carlo (MC) simulation is used to model the signal efficiency and the 
physics backgrounds. Simulation samples equivalent to approximately three times the accumulated data 
were used to model BB events, and samples equivalent to approximately 1.5 times the accumulated data 
were used to model continuum events where e + e~ — > uu, dd, ss, cc and t + t~. Three large samples of 
signal events are simulated, where a B + meson decays to e + v e (7.8 x 10 6 events), fi + u^ (7.6 x 10 6 events), 
or t + u t (19.4 x 10 6 events), and a B~ meson decays to an acceptable B mode. Beam related background 
and detector noise from data are overlaid on the simulated events. 

3 ANALYSIS METHOD 

Due to the presence of at least one neutrino in the final state, the B + — ► £ + V£ decay modes lack the 
kinematic constraints that are usually exploited in B decay searches to reject both continuum and BB 
backgrounds. The strategy adopted for this analysis is to reconstruct exclusively the decay of one of the B 
mesons in the event, referred to as the "tag" B. The remaining particle(s) in the event, referred to as the 
"signal B", are then compared with the signatures expected for B + — ► £ + U£. In order to avoid experimenter 
bias, the signal region in data is not examined ("blinded") until the final yield extraction is performed. 

The tag B is reconstructed in the set of semileptonic B decay modes B~ — > D°£~uX, where £ is 
e or fj, and X can be either nothing or a transition particle from a higher mass charm state decay, which 
we do not attempt to reconstruct (although those tags consistent with neutral B decays are vetoed). The 
D° is reconstructed in four decay modes: K~ir + , K~tt + it~it + , K~Tr + ir°, and K®ir + ir~ . The K® is 
reconstructed only in the mode K® — > ir + ir~. In a previous search for B + — > t + v t [10] we found that 
the low momentum transition daughter of D*° decays need not be reconstructed. Reconstructing the final 
state B — > D°£uX provides a higher efficiency but somewhat lower purity than the exclusive reconstruction 
method of B~ — ► D* Q £~~p£. In this analysis, we employ a technique, known as the "seeding" method, to 
recapture one photon from the X state. The decays B + — > e + v e and B + — > have not been previously 
searched for in the recoil of semileptonic tags. 

Since the r decays before reaching active detector elements, the B + — ► t + u t signal is searched for 
in both leptonic and hadronic r decay modes: r + — ► e + v e V T , r + — > /i + f M z7 r , r + — ► it + V T , and r + — > 
tt + tt°V t . The branching fractions of the above r decay modes are listed in Table 1. 
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Table 1: Branching fractions for the r decay modes used in the B + — > t + u t search [11]. 



Decay Mode 


Branching Fraction (%) 


t + — > e + ^ e I7 T 

r + — > 7r + F T 
r + — > 7r + 7r°F T 


17.84 ± 0.05 
17.36 ± 0.05 
10.90 ± 0.07 
25.50 ±0.10 



3.1 Tag B Reconstruction 

The tag B reconstruction proceeds as follows. First, we reconstruct the D° candidates in the aforementioned 
four decay modes using reconstructed tracks and photons where a tt° is included. The tracks are required 
to meet particle identification criteria consistent with the particle hypothesis, and are required to converge 
at a common vertex. The 7r° candidate is required to have invariant mass between 0.115-0.150 GeV/c 2 and 
its daughter photon candidates must have a minimum energy of 30 MeV. The mass of the reconstructed 
D° candidates in K~ir + , K~it + it~it + , and K®it + it~ modes are required to be within 20 MeV/c 2 of the 
nominal mass [11]. In the i\~~7r + 7r° decay mode, the mass is required to be within 35 MeV/c 2 of the 
nominal mass [11]; this wider mass window accounts for the BaBar detector' lower mass resolution when 
reconstructing particle candidates from neutral clusters, as opposed to reconstructing candidates involving 
charged tracks. 

Finally, D°£ candidates are reconstructed by combining the D° with an identified electron or muon with 
momentum above 0.8 GeV/c in the CM frame. The D° and £ candidates are required to meet at a common 
vertex. If more than one suitable D°£ candidate is reconstructed in an event, the best candidate is taken to 
be the one with the highest vertex probability. The uncorrected tag reconstruction efficiency in the signal 
MC simulation is 1.7% for B + -> t + v t , 1.1% for B + -> fi+u^ and 1.1% for B + -> e + v e . 

3.2 Selection of B + — ► £ + U£ signal candidates 

After the tag B reconstruction, in the signal B, we identify one of the following reconstructed particles: 
e + , fi + , 7r + , or p + . The e + and /i + can come from B + — > t + v t or directly from B + — ► or 
B + — > e + u e . Each signal B track must satisfy the following selection criteria: its point of closest approach 
to the interaction point is less than 2.5 cm along the beam axis and less than 1.5 cm transverse to the beam 
axis. 

The different reconstructed particles are assigned using a hierarchical selection involving kinematic 
constraints and particle identification. All of the signal decay modes for which we search contain only one 
track from the signal B. If more tracks that match the criteria stated above are present after a tag B has been 
reconstructed, the event is rejected. If the track from the signal B is identified as a kaon, it is rejected. Since 
we search for only one track in the signal B, we classify that track as one of the following in the priority 
given. 

• If the track satisfies particle identification as a muon, it is classified as such. 

• If the track satisfies particle identification as an electron, it is classified as such. We apply Bremsstrahlung 
radiation recovery techniques to identify as many electrons as possible. 
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• If the track can be combined with a 7r° to form a p + with a common vertex, it is classified as such. 
The invariant mass of a 7r° candidate must be between 0.115-0.150 GeV/c; the shower shape of the 
daughter photon candidates must be consistent with an electromagnetic shower shape, and the photons 
must have a minimum energy of 50 MeV in the CM frame. 

• If the track is not accepted by any of the above filters, it is classified as a ir + by default. 

Background consists primarily of B + B~ events in which the tag B meson has been correctly reconstructed 
and the recoil side contains one signal candidate track and additional particles which are not reconstructed by 
the tracking detectors or calorimeters. Typically these events contain K\ candidates and/or neutrinos, and 
frequently also additional charged or neutral particles which pass outside of the tracking and calorimeter 
acceptance. Background events also contain B°B° events. In addition some excess events in data, most 
likely from two-photon and QED processes which are not modeled in the MC simulation, are also seen. 

Multiple variables are used to suppress backgrounds. Most are combined into two likelihood ratios 
(LHRs), which are probability distributions designed to produce maximum separation between signal and 
background. Two variables are reserved for individual use due to their discriminating power. They are the 
momentum of a signal lepton in the B + rest frame (p sig e ) and the total energy recorded in the detector that 
is not assigned to the tag or signal B (E extTa ). 

Due to the presence of the neutrino in the products of the tag B, the direction of neither B can be 
known accurately. Instead, cos 9 B -D (> e ( tne cosine of the angle between the D°£ candidate and the B meson 
momenta) is calculated in the T(4S) rest frame. 

cos Q B -bH = 3 D i B DH , (2) 
2\PB\\PD°e\ 

where (E D o e , Vo^d and (Eb, Pb) are the four-momenta in the CM frame, and m D o e and mg are the masses 
of the D°£ candidate and tag B~ meson, respectively. Eb and the magnitude of ps are calculated from 
the beam energy: E B = -E'cm/2 and \ps\ = \J E 2 B — m? B . This definition assumes that the only missing 
particle in the tag B decay is a massless neutrino. Events in which the tag B daughters include a D and no 
higher mass charmed states are more common in the physical region, but other events have a larger tail into 
the non-physical region cos 9 B -d {> £ < — 1- 

For the reconstructed leptons in the signal B, we estimate the momentum of the signal lepton in the 
signal B + rest frame (p sig £ ) by averaging around the cone formed by cos 6 B -d°£- Since B + — > e + v e and 
B + — > n + Vn are two-body decays, for true signal events, p' sig e should exhibit a peak at 

If an event has a reconstructed signal muon candidate and p' sig e > 2.3 GeV/ c, it is classified as a B + — > 
fj, + u fl candidate; otherwise it is classified as a r + — > /i + z^F T candidate. If an event has a reconstructed 
signal electron candidate and p sig e > 2.25 GeV/c, it is classified as a B + — > e + u e candidate; otherwise it 
is classified as a r + — ► e + u e V T candidate. 

In an ideal B + — ► £ + U£ decay, we reconstruct all tracks and clusters associated with the real decay. The 
only unreconstructed particles would be neutrinos, which leave no energy in the detector. Therefore, we 
expect our signal to concentrate near zero .Bextra- We require a minimum energy of 30 MeV for any neutral 
cluster. 

After the D° has been reconstructed, a "seeding" algorithm adds a photon (called the "seed photon") to 
the reconstructed D° and reevaluates cos 6 B -d"£- The seeding algorithm performs this procedure with all 
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photons that do not overlap with the tag B and have CM energy less than 300 MeV. If a seed photon causes 
cos 9 B _ D o£ to become closer to (but not greater than) 1, it is selected. We seek to move cos9 B _ D o i closer 
to 1 because events containing real D* mesons usually appear in the low tail of the cos 9 B _ D oi distribution. 
If more than one photon satisfies these conditions, the one which moves AM = m D *o — m D o closest to the 
nominal value of 142.12 MeV/c 2 [11] is used. The photon is removed from E'extra and the event is stored 
with the modified -Eextra and cos 9 B _ D o i variables. 

We use a single photon to account for decays such as D*° — > D°j and D*° — > D°7r° (ir° — > 77). We 
studied the possibility of including a second seed photon, but it did not produce a significant improvement 
in performance. 



3.3 Likelihood Ratios 

To take advantage of shape differences between variables, we use two LHRs that consist of several probabil- 
ity density functions (PDFs). Separate LHRs are generated for BB and continuum background suppression. 
Two PDFs are generated for each variable in a LHR. One uses the signal MC sample and is treated as a prob- 
ability P s (x), where x is the value of the PDF variable. The other uses the relevant generic MC samples and 
is treated as a probability Pb{x). These two PDFs are combined to form a probability Pi(x) 

where each i represents a different variable. Bins in a distribution that are more likely to contain background 
events have Pi(x) closer to 0; bins that are more likely to contain signal events have Pi(x) closer to 1. Each 
LHR is formed by multiplying all P{(x) together: 

LHR(x) = n^W- (5) 

i 

Ideally, a LHR is a doubly peaked distribution with background events forming a peak near zero and 
signal events forming a peak near 1. The PDFs are created using MC samples with all tag selection criteria 
applied. Any given event will have one LHR for BB events and one for continuum events, where the PDFs 
are selected based on the reconstructed decay mode of that event. 



3.3.1 Variables Included in LHR 

Multiple variables were considered for inclusion in the LHRs. For each of the 14 LHRs (7 decay modes x 
2 background types), a fixed signal yield was chosen. Each LHR was tested using signal MC samples and 
the background MC samples it was designed to reject. The test were performed on MC after the appropriate 
decay mode was selected and with E extra required to be less than 1.5 GeV. For a given LHR, a baseline 
performance was calculated by using all prospective variables. A cut was placed on the LHR in question to 
produce the chosen signal yield, and the Punzi Figure of Merit (FOMp unz i) was calculated [12]. 

FOMpunai = N J;^ , (« 

where A s j g is the signal yield and Abg is the background yield. N a is the number of standard deviations 
desired from the result. We use N a = 3. 

The Punzi FOM is better suited to searches for small signals on small backgrounds. It is designed to 
prevent optimization algorithms from reducing the background to zero and creating a undesirably low signal 
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yield. We used it for the LHRs for all modes because significance (-/V sig / y (iV s ig + Nbg)) sometimes caused 
our optimization algorithm to reduce the background to zero. 

Each PDF was tested by removing it (and only it) from the LHR. The LHR was scanned again until the 
chosen signal yield was reached, and the FOM was recalculated. If removing the PDF increased the FOM 
by a statistically significant amount, it was not included in the final LHR. We also removed any variables 
that should not, for physical reasons, improve the analysis. For instance, variables related to the tag B would 
not reject BB background because the vast majority of it has a properly reconstructed tag B. The results of 
this selection process are shown in Table 2. Four of these variables were combined to form 2-D PDFs. The 
remainder were used as 1-D PDFs. 

• Separation Between the Signal and Tag B Vertices (Az/aAzY- 

Due to the neutrinos on both sides of the event, the vertices of the reconstructed signal tracks and 
neutral clusters do not correspond exactly to the true B decay points. However, the reconstructed 
vertices are still displaced in space while tracks from continuum processes tend to point back to the 
interaction point. 

We calculate the displacement between the putative B vertices divided by the uncertainty on that 
displacement. Continuum events are distributed more strongly towards zero than true BB events. 

• Net Event Charge: 

In our previous search [10], we noted a drop in tag efficiency over the lifetime of the experiment. We 
found that approximately half of this drop was due to a requirement that the event have zero net charge 
in order to pass tag selection. To avoid this drop in efficiency while retaining the discriminating power 
of this variable, the net charge of the event was moved from tag selection to the LHRs. 

• Ratio of the Second to Zeroth Fox- Wolfram Moment (i?2Aii): 

The Fox-Wolfram moments are rotationally invariant kinematic quantities designed to quantify the 
shapes of events resulting from e + e~ collisions. They are denoted Hi, where / is the number of the 
moment. 



Hi 



21 + . 

171= — / 



V s 



(7) 



where i runs over all hadrons in the event, pi are the momenta of the hadrons (in the CM frame), 
Y" z m are the spherical harmonics, 9 is the angle of the momentum with respect to the z axis, and y/s 
center-of-mass energy of the collision [13]. This variable is the ratio 

#2 A11 = ^. (8) 

We do not place a cut on this variable. We use it as a PDF in the LHRs of several modes, just as we 
use all of the other variables in this list. 

cos 9 B -D°e '■ 

This angular variable is defined in Section 3.2 and Equation 2. 
D° Decay Mode: 

Each D° decay mode is assigned an integer, and these integers form a distribution for an MC sample. 
Signal and background events have different distributions of D° decay mode. For instance, true signal 
events are found recoiling against aD ^ K~tt + more often than BB background events. 
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• Center of Mass Momentum for the Tag K (p^ a „ K -) and t (p£w-) 



• Tag K Selector: 

If the tag D° decay produces a putative charged K, we assign an integer value to the track correspond- 
ing to level of K particle identification passed by the track, with increasing values indicating tighter 
selection criteria. This is an integer value ranging from 10-14; each track is assigned to exactly one of 
these numbers. The value 16 is assigned to candidates. Signal events concentrate more strongly 
than background events at higher values. 

• Minimum Invariant Mass of Any 2 Reconstructed Tracks (Mf 1111 ): 

Since the minimum invariant mass of any three tracks M™" 1 was a useful variable in our previous 
search [10], we decide to try using M™ 11 . As the name suggests, it is the smallest invariant mass 
produced by any combination of two tracks used to reconstruct the signal B. 

• rn u : 

If an event contains two putative lepton tracks, their invariant mass is calculated and stored as mu- 
This variable was originally developed to remove pair-produced leptons. It is obviously highly corre- 
lated with Mf 1111 , so both variables are never used in the same LHR. 

• Signal n Selector: 

If the track from the signal B decay passes the particle identification requirements to become a puta- 
tive n + , we check if it passes a stricter level of fi + identification. If so, the event is assigned the value 
one for this PDF. Otherwise, it is assigned the value zero. More Continuum background than signal 
accumulates at zero. 

• Signal K + Selector: 

We wish to suppress the misreconstruction of K + as pions or leptons from the signal B, so we include 
K selection in the LHR. Specifically, this PDF is set to zero if a signal track passes loose K + particle 
identification requirements. 

• Vertex Status: 

For those r modes that involve neutral clusters (i.e r + — ► 7r + 7r F T ), a vertex is created. The quality of 
that vertex is reported as an integer from zero through four, which is included as PDF. Zero indicates 
that the vertex fit was successful. Other values indicate various failure modes for the fit. The only 
failure mode that occurs in this analysis is that the fit does not converge 

• Reconstructed Mass of the r Daughter (m p +): 

The decay r + — > 7r + 7r°z7 r often proceeds through the p + resonance. For true signal events, a peak 
at the resonance mass appears in the invariant mass distribution of the signal track and neutrals. 
Background events yield a flat or linear distribution. 

• Center of Mass Momentum for the ir + and ir° in r + — ► ir + ir°V T P*o)- 

2-D PDFs are two-dimensional histograms that contain two variable distributions. We use 2-D PDFs in 
cases where we want to exploit two variables that are highly correlated or that have a stronger separation 
when combined than when separate. 



14 



Table 2: This is a list of all variables used as signal PDFs for the various decay modes. Variables in bold are 
used for both Continuum and BB background. 



t + — > e + v e v T 




T + — » 7T + ^ T 






B+ -> e+!/ e 


Net Charge 

-R2 a n 

PTagf 

Tag K~ Sel. Level 
D° Dec. Mode 


COS 9 B _ D o e 






Az/lTAz 

m u 

^'T'unreco'dir 

CosrY-pY 


Az/(TAz 

_ min 

m 2 

f^unreco'dir 

Signal Sel. 
Signal fj, Sel. 
CosrY-pY 


CosrY-pY 


Az/(7Az 

Signal K+ Sel. 
Vertex Status 


Az/aAz 
m u 

^^unreco'dir 

Signal fj, Sel. 


^T'unreco'dir 



• "T-unrcco vs. cos(p^ lss ) (m unrcco -dir): 

The total invariant mass and initial momentum of each event are well known from beam information. 
Since neutrinos escape undetected from each event, we expect the total reconstructed invariant mass 
to be less than what the beam provides. The difference is called the unreconstructed mass (m unreC o)- 
The missing momentum (p™ lss ) is similarly defined, where 9 is the angle with respect to the beam 
line. 

The PEP-II beam pipe corresponds to values of cos(p™ lss ) near ±1. One source of background is 
events in which real particles are lost down the beam pipe, which is outside of the detector coverage. 
Since they are not reconstructed, they can be misinterpreted as neutrinos. This PDF allows the LHR to 
account for background events that have high unreconstructed mass but are likely to have lost particles 
down the beam pipe. 

• cos^_ y vs. \p' Y \ (CosrY-pY): 

cos 9' t _y is the equivalent of cos 9 B-D°e f° r tne signal B. Y represents all of the reconstructed 
daughters of the signal B, and 9' T _ Y is the calculated angle between Y and r momenta in the signal 
B rest frame. \p' Y | is the estimated momentum of Y in the signal B rest frame. It is the same as 
p sig t defined to include the hadronic r decay modes; it is calculated with the same average around 

the cone formed by cos 9 B _ D o£. Since p' sig e is such a powerful variable for selected B + — > e + u e and 
B + — ► we tested it for B + — > t + u t . We found that it was not very useful unless combined in 

this PDF. This PDF is not used in B + — ► e + v e or B + — > ^ + z^ reconstruction because p sig e is kept 
as a separate variable. 
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3.4 Optimization 



We use three (-E ex tra> ^HR^, and LHR cont ) variables in our final selection of the five r decay modes. For 
B + — > e + u e and B + — > /U+i^, we add a fourth variable p' sig e . The final requirements for these variables 
are obtained by optimizing on a figure of merit (FOM). For the r decay modes, we choose the FOM to be 
significance (N s \ g / (N sig + iVec))- F° r the other two leptonic B decay modes, we use the Punzi FOM 
(Equation 6). 

In the B + — > r + z/ r mode, both signal and background are large enough that optimizations perform well 
with standard significance. We make this determination based on the iV s i g expected from the SM predictions 
for the branching fractions. The optimized selection criteria are shown in Table 3. 

Table 3: Optimized ranges from which we accept signal candidates, which are mostly given by our opti- 
mization procedure. The exceptions are the p sig £ ranges for the two leptonic r decay modes, which were 
chosen to separate the leptonic r decays from B + — > e + u e and B + — > jj, + v^. 



Mode 


-E"extra 


LHR BB 


LHR cont . 


> 


r+ 




[0,0.24] GeV 


[0.74,1] 


[0.16,1] 


[0.00,2.25] GeV/c 


r+ - 




[0,0.24] GeV 


[0.14,1] 


[0.72,1] 


[0.00,2.30] GeV/c 


r+ 


TT + V T 


[0,0.35] GeV 


[0.57,1] 


[0.80,1] 




r+ - 


ir + ir°V T 


[0,0.24] GeV 


[0.97,1] 


[0.95,1] 




5+ -. 




[0,0.72] GeV 


[0.33,1] 


[0.75,1] 


[2.45,2.92] GeV/c 


B+ - 




[0,0.57] GeV 


[0.00,1] 


[0.01,1] 


[2.52,3.02] GeV/c 



After we examined the data in the signal region, we discovered an excess of data above our MC simula- 
tions at low values of mu, which is the minimum invariant mass of any two leptons. These events constitute 
an unmodeled background and are most likely due to photon pair conversion in the material of the detector. 
We decided remove all events below a certain value of mu after all other analysis cuts had been applied. 
This value was chosen using only signal and background MC simulations with the optimization technique 
described in this section. The result excludes events in the r + — > e + u e V T channel with mu < 0.29 GeV/ c 2 , 
which constitutes less than two percent of our signal MC sample. All efficiencies and yields in this note 
include the effects of this requirement. 

3.5 Signal Efficiency 

The signal B selection efficiencies for the decay modes are determined from signal MC simulation and 
summarized in Table 4. The signal efficiencies correspond to the number of events selected in a specific 
signal decay mode, given that a tag B has been reconstructed. 

The selection efficiency for r + — ► ^ + v^U T is low compared to that of the r + — > e + zv e !7 r mode because 
the momentum spectrum of the signal muons peaks below 1.2 GeV/c, where the muon detection efficiency 
is low. Since no minimum momentum requirement and no tight pion identification criteria are applied to the 
r + — > tt + V t signal selection, electron and muon signal tracks that fail particle identification requirement 
get selected in this mode. Any true r + — ► ir + ir°V T signal events, with a missed 7r° are also included in 
r + — > 7r + I7 T selection mode. Therefore, the r + — > tt + 17 t selection mode has the highest signal efficiency. 
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Table 4: Overall efficiency {e = e s i g x e tag ) of optimized signal selection for all modes before systematic 
corrections. 



Mode 


£ sig 


e(xl0" 4 ) 


t + — > e + f e I7 T 

T + — > 7T + F r 

r + — > 7r + 7r°F r 


(1.987 ±0.043)% 
(1.610 ±0.038)% 

(2.48 ±0.05)% 
(0.859 ± 0.028)% 


3.38 ± 0.07 
2.73 ± 0.06 
4.21 ± 0.08 
1.46 ±0.05 


S+ -► T+Z/ T 


(6.94 ±0.08)% 


11.78 ±0.13 




(30.92 ± 0.36)% 


32.54 ± 0.36 


S+ -► e+i/ e 


(36.98 ± 0.38)% 


40.43 ± 0.40 



3.6 Background Estimation from E extra Sidebands 

We define the "sideband" region as .Eextra > 0.6 GeV, except for B + — > where it is defined as 

^extra > 0.72 GeV. The "signal region" is defined separately for each signal mode using the optimized cuts 
on E'extra given in Table 3. Distributions of -Eextra for the signal decay modes are shown in Figures 3 - 6(b). 

For each mode, after applying the optimized final selections (except E'extra), the number of MC events 
in the signal region (-/V MCiSig ) and side band (A?MC,SideB) are counted and their ratio (R MC ) is obtained. 

^MC _ ^MC.Sig 
^MCSideB 

Using the number of data events in the side band (-/Vdata,sideB) and the ratio R MC , the number of expected 
background events in the signal region in data (A^x P ,sig) is estimated. 

-^exp.Sig = -Ndata,SideB ' R M ° 

Table 5 shows the background predictions from the -E e xtra sideband. We verify that the background 
predictions given by this sideband are consistent with the D° mass, LHR cont , LHR B -^, and p sig e side- 
bands, where applicable. We also studied the predictions given when the Eextra sideband is loosened to be 
> 0.8 GeV or tightened to include all events but the signal region in each mode. These predictions are also 
consistent with Table 5. 



4 SYSTEMATIC STUDIES 

The branching fraction for any of the decay modes in this analysis is given by 

B{B+ -> e+ue) = ^° bs ~ Nbg , (9) 

where N Q ^ S is the total number of events observed in the signal region, Nbg is the predicted number of 
events from background in the signal region; we use the values given as iVe X p,sig in Table 5 for Nbg- By 
definition = N s - lg + Nbg- N bb - * s tne tota ^ num b er of T(4S) decays in the data set, and the efficiencies 
can have different values for each mode. Each of these variables, except N ^ s , brings a systematic error into 
the branching fraction. 
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Table 5: BG Predictions from -E e xtra sideband. Rmc is the ratio of events in the sideband to events in the 
signal region of E'extra in the background MC. AT dataj g ideB is the number of events in the -E e xtra sideband in 
data. iVMC,Sig is the number of normalized events in the -Eextra signal region of the background MC samples. 
This is the background prediction taken solely from the MC samples. iV C xp,Sig is the product of Rmc and 
iVd a ta,SideB '■> it is the background prediction extrapolated from the data sideband using the MC samples. 



Mode 




^Vdata,SideB 


-^exp,Sig 


■^MC.Sig 


r + — > e + f e F T 


0.322 ± 0.040 


284 ± 17 


91 ± 13 


98 ± 11 


T+ -> fJ, + V,J7 T 


0.128 ±0.012 


1070 ± 33 


137 ± 13 


136 ± 12 


T + — ► 7T + F T 


0.033 ± 0.003 


6990 ± 80 


233 ± 19 


212 ± 17 


T + — > TT + TT V T 


0.035 ± 0.005 


1684 ± 41 


59 ±9 


62 ±9 


B+ 


1.1 ± 0.6 


14.0 ± 3.7 


15 ± 10 


12 ±5 


B+ -► e+z/ e 


0.57 ± 0.25 


42 ±6 


24 ± 11 


15 ±5 



4.1 Double Tag Control Sample 

To assess the agreement between Data and MC samples, we used two sets of control samples: sidebands 
and "double tagged" events. Events where both of the B mesons are reconstructed in tagging modes, 
B~ — > D(*)°£~V£ vs. B + — > D^°£ + V£, are referred to as "double tag" events. Due to both the large 
branching fraction of DiuX decays and the high tagging efficiency for reconstructing these events, a sizable 
sample (« 3.4 x 10 6 ) of such events are available in the on-resonance dataset. 

For double tag events, we first applied the tag selection requirements described in section 3.1 to both 
of the tag candidates. This procedure resulted in noticeable shape and yield discrepancies between data 
and MC, as seen in Figure 2(a). In order to improve agreement, several additional selection criteria were 
imposed. These are based on the selection criteria from the previous search for B + — ► t + v t [10]. We 
require that —2.0 < cos9 B _ D o£ < 1.1 for both B decays and the event has zero net charge. The resulting 
distribution is shown in Figure 2(b). The yield disagreement has improved, and the agreement in shape is 
better. To test how much of the disagreement is due to yield and how much is due to shape, we normalize 
data and MC to unit area with the same cuts as in Figure 2(b). The resulting distribution is shown in Figure 
2(c), which shows excellent agreement in the shape of the distribution. 

4.2 Systematic Error from Background Prediction (N B g) 

We use the ratio of data to MC samples in our background prediction, and the statistical error on that 
prediction is already large. The shape of our data and MC samples agree well in the sideband region, so we 
have no need to apply an additional systematic error. Therefore, we accept this as the total error and apply 
no further systematic correction or uncertainty. 

4.3 Systematic Error from B-Counting (N BB ) 

The estimation of the number of B + B~ events present in our data sample has a small uncertainty of 1.1% 
[14]. 
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Figure 2: Distribution of E sxtISl in double tags with different sets of selection requirements. The Data and 
MC samples are normalized to Data luminosity in 2(a) and 2(b). In 2(b), additional analysis requirements 
are applied in order to improve agreement. In 2(c) the same requirements are applied as in 2(b), but Data 
and MC are normalized to unit area. The gray rectangles represent the extent of the error bars on the MC 
histogram. 
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4.4 Systematic Error from Tagging Efficiency (£t ag ) 

We pursue a procedure that attempts to combine what we know about both the single- and double-tagged 
samples. We define the efficiency 

e 2 = N 2 /N 1 , 

where Ni is the single-tagged events, and N 2 is the number of double-tagged events, as defined in Section 
4.1. We use the ratio of e 2 in Data and our MC samples as a systematic correction to the tag efficiency 
(0.891). We take the uncertainty on this ratio as a systematic uncertainty (0.021). These numbers were 
extracted using only those events reconstructed as D° — > K~ir + for the first tag. We verify that this 
correction is consistent with the correction calculated using D° — ► K~it + it~it + on the first tag. 

4.5 Systematic Error from Signal Efficiency (e sig ) 

The systematic uncertainties on the signal selection efficiency for each signal mode have contributions from 
©extra mis-modeling, tracking efficiency, particle identification, and 7r° reconstruction for the r + — > 7r + 7r°I7 r 
mode. 

The systematic correction and uncertainty on the © ex tra shape is taken from the double tag © ex tra distri- 
bution described in Section 4.1. The double tagged events provide us with a means of comparing data and 
simulation, using an independent control sample, to extract this uncertainty. We use the © ex tra distributions 
shown in Figure 2(b) to extract the yield of candidates satisfying -©extra < 0.6 GeV. This yield is then com- 
pared to the number of candidates in the full distribution. Comparing the ratio extracted from MC to that 
extracted from data yields a correction factor, the error on which is taken as the systematic uncertainty for 
©extra- We extract a correction and uncertainty of 1.015 ± 0.021. 

Since the particle identification algorithms have not changed since our previous search for B + — ► t + u t 
[10], we use the same values as conservative estimates of our current systematic errors. Since the previous 
analysis only used the four single-prong r decay modes, we must extrapolate. For the r + — ► tt + -k~it + V t 
mode, we apply the -k correction three times and take triple the uncertainty as our a\ uncertainty. For the 
B + — > /ti+z^ and B + — ► e + v e modes, we use the [i and e systematics, which is a conservative choice 
because particle identification is more effective at the high momenta that characterize two-body B decays. 

All multiplicative contributions to the systematic uncertainty are summarized in Table 6. The corrected 
efficiencies are shown in Table 7. We extract a total multiplicative systematic uncertainty of 3.6% for 
B+ -> T +v T , 4.4% for B + -> fi+v^, and 4.0% for B + -> e + v e . 

5 RESULTS 

After finalizing the signal selection criteria, the signal region in the on-resonance data is examined. Table 8 
lists the number of observed events in on-resonance data in the signal region, together with the expected 
number of background events in the signal region. Figures 3 and 4 show the ©extra distribution in data and 
simulation for each of the r decay modes considered. Data is overlaid on the summed MC contribution, 
scaled to the dataset luminosity, and signal MC is plotted for comparison. Figure 5 shows the © cx tra distri- 
bution for all r decay modes combined, with MC scaled to the sideband data yield. Figure 6(a) shows the 
©extra distribution for B + — > with and without MC scaling to the sideband data yield. Figure 6(b) 

shows the © cx tra distribution for B + — > e + u e with and without MC scaling to the sideband data yield. 

We use the method developed by Feldman and Cousins [15], which is designed to produce an upper 
limit for null results and a two-sided confidence interval for non-null results. The Feldman Cousins method 
begins with the construction of a confidence belt, which is a two-dimensional histogram A tmc vs. A s ; g . 
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Table 6: Summary of Systematic corrections, uncertainties, and fractional uncertainties. Note that the 
numbers from the rightmost column of Table 7 are included in the rightmost column of this Table. 



Source 


Applicable Mode(s) 


Correction 


Fractional Uncertainty (%) 


B Counting 


All 


1.0 


1.1 


Tag efficiency 


All 


0.891 ± 0.021 


2.4 


-E'extra 


All 


1.015 ±0.021 


2.1 


7T° Reconstruction 


T + — > 7t + Tr V T 


0.984 ± 0.030 


3.0 


Tracking Efficiency 


t + — ► e + v e 17 T 


1.0 


0.36 




T+ — > /i + ty7 T 


1.0 


0.36 




T + — > 7T + I7 T 


1.0 


0.36 




r + — > 7r + 7r F r 


1.0 


0.36 




-► /i+l/u 


1.0 


0.36 




5+ -► e+zv e 


1.0 


0.36 


Particle Identification 


r + — ► e + ^ e I7 T 


1.01 


2.5 




r+ -> /i + i/ M F T 


0.92 


3.1 




r + — ► 7r + F T 


1.02 


0.8 




T + — > 7t + TT°V T 


1.00 


1.5 




B+ 


0.92 


3.1 




5+ -► e+^ e 


1.01 


2.5 



A^ s ig represents the extracted signal yield for an ensemble of experiments for each value of iV true . iV true is 
the actual number of signal events used as the central value to generate each ensemble of experiments. We 
generate these distributions using a random number generator. 

For each value of AT true , we generate two sets of 100,000 random numbers. For the first set, a Poisson 
random number generator is used with the central value set to iVt rue . The second set is based on the back- 
ground predictions from Table 5. A Gaussian random number generator is used with the center set to the 
central value from the table (Nbg) an d the width equal to \/Nbg- This value is then used as the center of 
another random number distribution (iVbackX with a width equal to the error on (Nbg) from Table 5. This 
procedure is used to account for both the error on the background prediction and the statistical error on the 
total number of observed events (iV bs)- These two numbers are summed for each experiment to form iV bg. 

A^sig = A^obs - AW (10) 

To smooth the statistical fluctuations resulting from the random number generation, we fit the distribu- 
tion for each value of A truc to the sum of two Gaussians. Asig is then reassigned its value from the fitting 
functions. For each A true , we must define an acceptance region in A^ig that will determine our upper limit 
or central value with uncertainty. The Feldman Cousins method defines this acceptance region without re- 
ferring to data or any bias regarding whether we seek an upper limit or branching fraction. For each bin, we 
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Table 7: The corrected tag and signal efficiencies. Two errors are quoted: the first is the MC statistical 
uncertainty, and the second is the systematic error computed from the sources in this section. 



Efficiency 


Corrected 


Fractional Systematic Error (%) 


Tag(£+ -► t + u t ) 


(1.514 ±0.003 ±0.036)% 


2.4 


Tag(B+ n+vj 


(0.937 ±0.003 ±0.022)% 


2.4 


Tag(5+ -► e+i/e) 


(0.974 ±0.003 ±0.023)% 


2.4 


(T + -»e+I/ e F T ) 
sig 


(2.04 ±0.04 ±0.07)% 


3.3 


(t+-^+iv7 t ) 
sig 


(1.50 ±0.04 ±0.06)% 


3.7 


(r+^7r+j/ T ) 
sig 


(2.57 ±0.05 ±0.06)% 


2.2 


£ sig 


(0.86 ±0.03 ±0.03)% 


4.0 


£ sig 


(6.97 ±0.08 ±0.20)% 


2.8 


sig 


(28.9 ±0.3 ±1.1)% 


3.7 


sig 


(37.9 ±0.4 ±1.2)% 


3.3 



calculated the ratio 

_ P(jV sig |iVt rU e) 

~P(iV sig |iV best )' ; 

where iVbcst is the value of N truc that maximizes the probability of observing iV s ig. Thus, R ranges between 
and 1. For each value of AT true , we sort the bins of N sig in order of descending R. The probabilities 
P(iV s ig|JVt rue ) are summed in this order until the desired confidence level is reached. The resulting distribu- 
tions are shown in Figure 7. We calculate iV s ig and draw a vertical line through the appropriate distribution 
at iV s ig. The upper and lower limits are determined by the intersection of that line with the appropriate 
confidence bands. 

We determine the following preliminary branching fraction: 

B{B + t + v t ) = (1.8 ± 0.8 ± 0.1) x 10" 4 . (12) 

Using the Feldman-Cousins method, we determine this result excludes the null hypothesis at the 2.4cr level, 
including all systematic uncertainties. 

We notice that the electron mode has a much larger excess than the other r decay channels. We per- 
formed several tests and cross-checks to determine if this discrepancy is a statistical fluctuation or due to 
a background mode not included in our MC simulation. We have tested for the following potential back- 
ground contributions: two photon fusion QED events, "events" that contain two overlapping e + e _ colli- 
sions, overzealous Bremsstrahlung recovery, photon pair production where the e + and e~ are reconstructed 
as the tag and signal lepton, and photon pair production events where one lepton is lost and the other is 
reconstructed as the signal electron. 

In two-photon fusion events, the e + and e~ each emit a photon; the two photons interact to to produce 
multiple hadrons. When these hadrons are misreconstructed as a D° and the e + and e~ are identified as a 
tag and signal lepton, the result is a signal-like event with low extra energy. Since the signal lepton is an e^, 
these events will only populate the electron r decay channel. Two-photon fusion events rarely produce real 
D° mesons, so they should populate the sidebands and peak of the D° mass distribution. However, we find 
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Figure 3: Total extra energy is plotted after all cuts have been applied in the mode (a) r + — > e + v e V T and 
(b) t + — > ^ + u^U T . The background MC samples have been scaled according to the ratio of predicted 
backgrounds from data and MC as presented in section 3.6 and summed together. The grey rectangles 
represent the extent of the error bars on the MC histogram. The signal region consists of the region of E extm 
to the left of the vertical dashed line. Simulated B + — > t + u t signal MC is plotted (lower) for comparison. 

that the background predictions taken from he D° mass sidebands (123.8 ± 17.7 events) is consistent with 
the prediction from the -Eextra sideband (Table 8). In the -Eextra distribution for r + — > e + v e V T taken from 
the sidebands of the D° mass distribution, the data agree will with the MC simulation. This shows that the 
excess of signal events in this mode occurs mostly in the peak of the D° mass distribution. Since two-photon 
fusion events rarely produce real D°, any excess from this process should also be in the sidebands. No such 
excess was present, so we can rule out this potential background. 

We examined the distribution of Az, which is the separation between the putative B vertices. Overlap- 
ping events should have wider separation than real events, and thus we should see an excess of data above 
the MC at high values for Az. No such excess was found. Another suggested source of background in 
this channel is Bremsstrahlung recovery that assigns more photons to an electron than it actually generated. 
This could move events into the signal range of E extm undeservedly. To test for this, we compare the E extm 
distributions for electrons with and without Bremsstrahlung recovery. When we turn off Bremsstrahlung 
recovery, we do so for electrons in both the tag and signal B. The comparison showed that data and MC 
simulation have very similar shapes regardless of whether Bremsstrahlung recovery is used. No suspicious 
excess appears when Bremsstrahlung recovery is activated, so this is not a likely source of the excess. 

To remove events in which both daughters of a photon pair conversion are reconstructed as the tag and 
signal leptons, we require that mil > 0.29 GeV/c 2 in the r + — > e + u e V T channel, as described in Section 
3.4. If a photon produces an e + e~ pair in the detector material, one of member of the pair is lost (e.g. down 
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Figure 4: Total extra energy is plotted after all cuts have been applied in the mode (a) r + — ► 7r + z7 T and (b) 
r + — > 7T + 7r°z7 T . The background MC samples have been scaled according to the ratio of predicted back- 
grounds from data and MC as presented in section 3.6 and summed together. The grey rectangles represent 
the extent of the error bars on the MC histogram. The signal region consists of the region of i^extra to the left 
of the vertical dashed line. Simulated B + — > t + v t signal MC is plotted (lower) for comparison. 

the beam pipe), a D ^ candidate is reconstructed, and these candidates pass all analysis cuts, we would 
reconstruct this as a false B + — » t + v t decay. However, either member of the e + e~ pair is equally likely 
to be lost, so any excess from this source of background would appear in events where the tag and signal 
lepton have the same electrical charge. However, the MC simulation matches the data for same-sign events 
very well. 

To test the probability of a statistical fluctuation producing the excess seen in the r + — > e + u e T' T channel, 
we fill a histogram with the branching fraction calculated from each r decay channel separately. The values 
and uncertainties are taken from the Feldman-Cousins method; systematic uncertainties are not included. 
We fit a constant to the branching fractions; the result is (1.8 ± 0.6) x 10~ 4 . The fit has a reduced x 2 of 
1.64, which corresponds to a probability of 18%. The results of this study are shown graphically in Figure 
8. 

Using \ V u b\ = (4.43 ± 0.54) x 10 -3 , we extract a preliminary value fs of 

f B = 230±57MeV. (13) 

The BABAR Collaboration previously published a statistically independent measurement of the B + — > 
t + v t branching fraction using tag B mesons decaying into fully hadronic final states[16]. We measured 

B(B + -» t + u t ) = (X.8±J;g) x 10~ 4 . (14) 
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Figure 5: Total extra energy is plotted after all cuts have been applied with all B + — > t + v t modes com- 
bined. Events in this distribution are required to pass all selection criteria. In addition, the background MC 
samples have been scaled according to the ratio of predicted backgrounds from data and MC as presented 
in section 3.6 and summed together. The grey rectangles represent the extent of the error bars on the MC 
histogram. Simulated B + — > t + v t signal MC is plotted (lower) for comparison. 

Combining these two measurements using a simple error-weighted averaged yields 

B(B + -► t + v t ) = (1.8 ± 0.6) x 1(T 4 , (15) 

which excludes zero at the 3.2<t level. 

Since the result for B{B + — > t + v t ) is still consistent with zero, we set an upper limit for all three 
modes at the 90% confidence level. 

B(B + -► t + u t ) < 3.2 x 10" 4 , (16) 

B(B + -► < 11 x 10~ 6 , (17) 

B(B + e + u e ) < 7.7 x 10~ 6 . (18) 

The upper limits for B(B + — > and B(B + — > e + i/ e ) are consistent with previous measurements 

including the current PDG values of < 1.7 x 10~ 6 and < 9.8 x 10~ 7 , respectively. The BaB/R Collaboration 
has set the latest limit on B{B + — > < 1.3 x 10 -6 [17]. The limits reported in this analysis are much 

higher than the best available limits because the tagging method produces a very low background at the cost 
of a low efficiency. Low backgrounds are more conducive to discovery; however, until we have enough data 
to make a statistically significant discovery, other methods will produce stricter upper limits. 
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Figure 6: Total extra energy is plotted after all cuts have been applied to B + — > /i + f M (left) and 1? + — > 
e + f e (right). Events in this distribution are required to pass all selection criteria. The background MC 
samples have been scaled according to the ratio of predicted backgrounds from data and MC as presented 
in section 3.6 and summed together. The grey rectangles represent the extent of the error bars on the MC 
histogram. The signal region consists of the region of E e ^ m to the left of the vertical dashed line. Simulated 
signal MC is plotted (lower) for comparison. 



6 CONCLUSIONS 

We have performed a search for the decay process B + — > £ + V£. To accomplish this, a sample of semilep- 
tonic B decays (B~ — > D^^l^vX) has been used to reconstruct one of the B mesons and the remaining 
information in the event is searched for evidence of B + — > 

For B + — > t + v t , we measure a preliminary branching fraction of 

B(B + -> t + u t ) = (1.8 ± 0.8 ± 0.1) x 10~ 4 , 

which excludes zero at the 2Aa level. 

We set preliminary upper limits at the 90% confidence level of 

B{B + -> t + u t ) < 3.2 x 10~ 4 , 

B(B + -» v+v^) < 11 x 10" 6 , 
B(B + -» e + v e ) < 7.7 x 10" 6 . 
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Figure 7: The confidence band produced by the Feldman Cousins method. The central band represents the 
la confidence level; the next band out represents the 90% confidence level. The vertical red line represents 
the value of iV s ig based on the number of observed events. Figure (a) is for B + — ► t + u t ; (b) is for 
B + — ► fi + u fl , and (c) is for B + — > e + u e . 

7 ACKNOWLEDGMENTS 

We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the excellent lu- 
minosity and machine conditions that have made this work possible. The success of this project also relies 
critically on the expertise and dedication of the computing organizations that support BaBar. The collab- 
orating institutions wish to thank SLAC for its support and the kind hospitality extended to them. This 
work is supported by the US Department of Energy and National Science Foundation, the Natural Sci- 
ences and Engineering Research Council (Canada), the Commissariat a l'Energie Atomique and Institut 
National de Physique Nucleaire et de Physique des Particules (France), the Bundesministerium fur Bildung 
und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionale di Fisica Nucle- 
are (Italy), the Foundation for Fundamental Research on Matter (The Netherlands), the Research Council 
of Norway, the Ministry of Education and Science of the Russian Federation, Ministerio de Educacidn y 
Ciencia (Spain), and the Science and Technology Facilities Council (United Kingdom). Individuals have 
received support from the Marie-Curie IEF program (European Union) and the A. P. Sloan Foundation. 



27 



Table 8: The observed number of on-resonance Data events (iV obs ) in the signal region is shown, together 
with the number of expected background events (-/Vbg), the corrected overall selection efficiency (e), and 
the branching fraction calculated from each channel. For the four t decay channels, the branching fraction 
shown is the branching fraction of B + — > r + v T using only this channel. All systematic uncertainties are 
included. 



Mode 


Expected 


Observed 


Overall 


Branching 




Background 


Events 


Efficiency (e) 


Fraction 




(JVbg) 


(N ohs ) 






t + — * e + v e V T 


91 ± 13 


148 


(3.08 ±0.14) x 10" 4 


(4.0 ± 1.2) x 10" 4 


r+ -> /x+i/ M F r 


137 ± 13 


148 


(2.28 ±0.11) x 10~ 4 


(l.0±l >) x 10- 4 


r + — * ir + V T 


233 ± 19 


243 


(3.89 ±0.15) x 10" 4 


o.e^i x io- 4 


T + — > 7T + 7r F r 


59 ±9 


71 


(1.30 ±0.07) x 10~ 4 


(2.0±{; 4 ) x 10" 4 




521 ± 31 


610 


(10.54 ±0.41) x 10~ 4 


(1.8 ±0.8 ±0.1) x 10~ 4 


B+ -» 


15 ± 10 


11 


(27.1 ± 1.2) x 10" 4 


< 11 x 10~ e @ 90% CL 


5+ -» e+z^ e 


24 ± 11 


17 


(36.9 ± 1.5) x 10" 4 


< 7.7 x 10~ s @ 90% CL 
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Figure 8: The branching fractions calculated from each r decay channel are plotted. 1 = r + — > e + v e v T , 
2 = t + — > ii+u^Vr, 3 = r + — ► 7r + I7 T , 4 = t + — > 7r + 7r I7 T . The horizontal line is a fit to a constant. 
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